arXiv:astro-ph/0503194v2 9 Mar 2005 


22nd TEXAS SYMPOSIUM on Relativistic Astrophysics, Stanford Univ., Dec. 13-17, 2004 


1 


The very early optical afterglow of GRBs, revealing the nature of the ejecta 


Y. Z. Fan 

Department of Physics, University of Nevada, Las Vegas NV 89154, USA 
Purple Mountain Observatory, Chinese Academy of Seience, Nanjing 210008, China, and 
National Astronomieal Observatories, Chinese Academy of Sciences, Beijing, 100012, China. 

Bing Zhang 

Department of Physics, University of Nevada, Las Vegas NV 89154, USA 

D. M. Wei 

Purple Mountain Observatory, Chinese Academy of Seience, Nanjing 210008, China, and 
National Astronomieal Observatories, Chinese Academy of Sciences, Beijing, 100012, China. 

We show that if the GRB ejecta itself is magnetized or neutron-rich, the very early afterglow of 
GRBs is very different from that powered by a pure ion fireball. In the Swift era, with the well 
monitored early afterglow data, we can potentially diagnose the ejecta composition and reveal the 
nature of the central GRB engine. 


I. INTRODUCTION 

Since the discovery of the optical flash of GRB 
990123 Q, great attention has been paid to the very 
earl y optic al afterglow (e.g., for theory see flTl 

f J3 EM EM EM 1^ 1^ ; for observation, see 
h4. 1^ 128 1 1 of soft, long GRBs. The 

very early optical afterglow of short GRBs and X- 
ray flashes have also been investigated in some detail 

iini. 

In nearly all these theoretical papers, the fireballs 
are assumed to be non-magnetized and neutron-free. 
However, in principle, the ejecta may be significantly 
magnetized (see |23l| and the references therein) or 
contain a large amout of neutrons (e.g. 0, @l)i or 

both (e.g. 113 )■ Recently we have studied the early 
afterglow emission (especially the reverse shock emis- 
sion) of a magnetized ejecta and a neutron-rich ejecta 
|9l ll2Ll32| . Here we summarize the main results. Com¬ 
paring these results with the early afterglow observa¬ 
tions, we may get insight into the initial composition 
of the outflow and reveal the nature of the central 
engine. 

II. REVERSE SHOCK EMISSION FROM 
MAGNETIZATED EJECTA 

A. Magnetized reverse shock emission from 
GRB 990123 and GRB 021211 

A bright optical flash was detected during the bright 
GRB 990123 (e.g., Q)- The peak V-band flux was 
8.9th magnitude. After the peak, the flux drops as 
^obs (before the peak, the flux increases as t^\^s). Such 
a sharp decrease has been modeled by the emission 
from the shocked electrons contained in the reverse 
shock region with adiabatic cooling |^. With 
a more careful investigation. Fan et al. found that 


by taking the reverse shock model and the physical 
parameters found in modeling the multi-wavelength 
afterglow from the forward shock emission, the theo¬ 
retical peak flux of the optical flash accounts for only 
3 X 10“"* of the observed value (~ IJy) Q- In order to 
remove this discrepancy, they suggested that the elec¬ 
tron and magnetic equipartition parameters, Ce and 
cb, should be 0.61 and 0.39, respectively. These are 
much different from the corresponding values for the 
late afterglow (~ 0.1 and ~ 0.001). 

In a more general discussion, Zhang, Kobayashi & 
Meszaros introduced a parameter TZb (which is de¬ 
fined as (cb/cb)^/^; where the superscripts r and / 
represent the reverse shock region and forward shock 
region, respectively) to trace the magnetization of the 
GRB outflows [33 • They found T^b ~15 for GRB 
990123 and 7^B > 1 for GRB 021211. 

These results have been also confirmed in several 
more detailed analysis of both bursts lililEl. As 
a result, it is quite robust to say that the reverse shock 
region is magnetized, at least for GRB 990123 and 
GRB 021211. In the following, we introduce a param¬ 
eter (T, the ratio between the electromagnetic energy 
flux and the particle energy flux, to describe the ini¬ 
tial magnetization of the outflow (we assume that the 
magnetic field is ordered). 


B. Reverse shock emission with mild 
magnetization 

As the outflow interacts with the surrounding 
medium, two shocks form. One is the forward shock 
expanding into the medium, and the other is the re¬ 
verse shock penetrating into the outflow. The forward 
shock jump condition is the same as that of a pure hy- 
drodynamical fireball model 0 , but the reverse shock 
jump condition is different (see equations (2-5) of |T3 
for the general form). 
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The novel features include: 1. At the reverse shock 
front, for cr > 0.01, the amplified ordered magnetic 
field (i^ord) is signihcantly stronger than that of the 
random one (i?ran)- In the ISM medium, the typical 
synchrotron radiation ly^ oc iJord is much lower than 
v-R. ~ 4.6 X lO^^Hz and the reverse shock-accelerated 
electrons are in the slow cooling regime. The observed 
flux scales as 

F,, oc oc (1) 

Therefore, the R-band reverse shock emission is 
stronger for a higher Rord (see FigCJ, which matches 
the observation of GRB 990123 and GRB 021211. 
However, for a > 0.1, the reverse shock is suppressed 
(in both and the ideal MHD approximation 
is taken. If magntic dissipation M is taken into ac¬ 
count, the result would be different), so that the re¬ 
verse shock peak level starts to decrease (see Fig^). 
In the wind medium case, the R-band reverse shock 
emission with magnetization is weaker (see Fig. EJ 
since the reverse shock is relativistic and the electrons 
are in the fast cooling phase, and r'm is above vr. The 
observed flux scales as 

F,.ocHordi^y'oci?-y", (2) 

where Vc oc B~^ is the cooling frequency. As a result, 
Fy decreaes with increasing a. 

2. In the very early afterglow phase, the outflow 
is ultrarelativistic. Due to the beaming effect, the 
area we view is very narrow where the orientation of 
magnetic field is nearly the same. The local high lin¬ 
ear polarization can not be averaged effectively. As 
a result, net high linear polarization is expected. By 
introducing a parameter b = Bo^d/B^an, the net linear 
polarization can be approximated by (e.g., |^) 

nnet«0.60^. (3) 

Even for a ^ 0.01 — 0.1, the linear polarization as high 
as 30% is expected. 


C. The reverse shock emission with arbitrary 
magnetization 

Zhang & Kobayashi have performed a detailed an¬ 
alytical investigation on the reverse shock emission 
with arbitrary magnetization |32|. For cr < 1, i.e., 
the mildly magnetized regime, their results are rather 
similar to that of Q- For cr > 1, i.e., the high-cr 
regime, it is found out that the reverse shock peak is 
broadened, mainly due to the separation of the shock 
crossing radius and the deceleration radius of the out¬ 
flow (see FigEJ. This novel feature can be regarded 
as a signature of high cr. 



FIG. 1: The very early R-band (rr = 4.6 x lO^'^Hz) light 
curve powered by the mildly magnetized outflow (the de¬ 
grees of the magnetization have been marked in the figure) 
interacting with the interstellar medium. The parame¬ 
ters taken in the calculation are: 2 = 1 , Akin = 10 ®®ergs, 
L = 2 X lO^^ergs s“^, p = 2.2, p = 300, n = lcm“®, 
£e = 0.3 and the radiation efficiency e = Ce, where Akin is 
the total initial energy of the outflow, L is the luminosity 
of the 7 —ray emission, p is the initial bulk Lorentz factor 
of the outflow, p is the power-law distribution index of the 
shock-accelerated electrons, and n is the number density 
of the medium. For a = 0 and the forward shock, it is 
assumed that cb = 0.01. From ref [^. 

Another important result obtained is that the sup¬ 
pression factor of the reverse shock in the strong mag¬ 
netic field regime is only mild as long as the shock 
is relativistic, and it saturates in the high-cr regime 
(see FigEj). This indicates that strong relativistic 
shocks still exist in the high-cr limit, which can ef¬ 
fectively convert kinetic energy into heat. The overall 
efficiency of converting jet energy into heat, however, 
decreases with increasing cr, mainly because the frac¬ 
tion of the kinetic energy in the total energy decreases 
1^. These results have been confirmed by Fan Wei 
& Zhang both analytically and numerically [2, see 
their equations (17-18) and figures 1 and 2. 

III. VERY EARLY OPTICAL AFTERGLOW 

LIGHTCURVES OF NEUTRON-FED GRBS 

In the rest frame of the ejecta, the free neutron has a 
mean lifetime ~ 900s. The corresponding /3—decay ra¬ 
dius reads Rp ^ 8x 10^^rn^2.5cm where Fn ^ 300 
is the bulk Lorentz factor of the neutrons (below we 
call it as N-ejecta). In the internal shock phase, the 
ion-ejecta (I-ejecta) has been decelerated, but the neu¬ 
trons are not. They move freely into the medium and 
decay into protons, neutrinos and electrons. These da- 
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FIG. 2: The very early R-band light curve powered by 
the mildly magnetized outflow (the degree of the magne¬ 
tization has been marked in the figure) interacting with 
a stellar wind. The parameters taken here are the same 
to those of figure 1 except n = 3 x 10^®i?“^cm“^. From 
ref 0 . 



FIG. 3: Sample early afterglow lightcurves for GRBs with 
an arbitrary magnetization parameter a. Following pa¬ 
rameters are adopted. iSkin = 10®^ergs, T = 20s, rj — 150, 
n = lcm~®, Eg = Ce = 0 . 1 , eb — 0 . 001 , p = 2 . 2 , and z = 1, 
where T is the duration of the GRB corrected by redshift. 
Both the forward shock and the reverse shock emission 
components are calculated and they are superposed to get 
the final lightcurve. Li ght curves are calculated for differ¬ 
ent a values. From ref 133 . 


cay products share their energy and momentum with 
the medium and form a mixture (the trail) moving 
with a bulk Lorentz factor about tens, the actual ve¬ 
locity mainly depends on the density of the medium. 

Comparing with that of the pure fireball or the 
magnetized fireball, the early afterglow of neutron-fed 



FIG. 4: The variations of six parameters, i.e., U 2 s, 7 is/ 72 i, 
e^ln^mpcf, n^jnx, P 6 . 2 /P 2 , and F, as a function of cr. 
The thick solid line is the Kennel-Coroniti solution m, 
denoting a 721 ^ 1 regime. Here U 2 s is the radial four 
velocity of region 2 measured in the shock frame; 71 s 
is the bulk Lorentz factoe of region 1 measured in the 
shock frame; 721 is the Lorentz factor of region 1 rela¬ 
tive to region 2 ; 62 is the thermal energy density of re¬ 
gion 2 (measured in the comoving frame); n 2 (ni) is the 
number density of region 2 ( 1 ) measured in its comoving 
frame; pb ,2 (P 2 ) is the magnetic field (thermal) pressure 
of region 2 measured in its comoving frame. The dashed 
lines, starting from the one closest to the thick line, are 
for 721 = 1000,100, 10, 5, 3, 1.5, respectively. Again the pa¬ 
rameters e 2 ln 2 vap(? and n 2 /ni are normalized to (721 — 1 ) 
and (4721 -I- 3), respectively. From ref liill . 


GRBs is very complicated. The detailed discussion 
has been presented in M Here we only summarize 
the main results. 

If the medium is a pre-stellar wind, the neutron 
trail moves slowly, mainly because the medium inertia 
is too large. The trail and the I-ejecta do not sepa¬ 
rate from each other, and a forward shock propagates 
into the trail directly. Three components contribute to 
the final emission, i.e. the forward shock, the reverse 
shock propagates into the Tejecta, and the unshocked 
trail emission. The latter is significant when y, the 
ratio of neutrons to protons, is large, since the inter¬ 
nal energy of the unshocked trail is large when the 
medium density is high. A typical neutron-rich wind- 
interaction lightcurve is a characterized by a promi¬ 
nent early plateau lasting for ~ 100 s followed by a 
normal power-law decay (FigEJ. We also show that in 
the wind case, the IC cooling effect due to the over- 
laping of the initial prompt 7 —ray with the shocks 
and the trail suppresses the very early R-band after¬ 
glow significantly. The neutron-fed signature is also 
dimmed (see Figl^b) and Fig[3(c) for a comparison). 

If the medium is a constant density ISM, part of the 
neutron decay products fall onto the medium, and the 
trail moves faster than the I-ejecta (In da , it is as- 
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FIG. 5: The early optical afterglow lightcurves of a 
neutron-fed long GRB in the wind interaction case, (a) 
The dynamical evolution of the LF of the shocked region 
as a function of time, (b) R-band lightcurves, with the 
IG cooling effect due to the prompt 7 —rays interacting 
with the shocked regions being ignored. Thick lines in¬ 
clude contributions from all emission components, includ¬ 
ing the FS, RS and the neutron decay trail. Thin lines are 
for trail emission only. The dotted, dash-dotted, dashed 
and solid lines represent x = 0-0, 0.1, 0.5, 1.0 respec¬ 
tively. Following input parameters are adopted in the cal¬ 
culations: Etot = 2.0 X lO^^ergs, A = lO^^cm, z = 1 
[i.e. dn = 2.2 x 10 ^®cm, which corresponds to the stan¬ 
dard {Urn, Ha) = (0.3,0.7) ACDM cosmological model], 
Fn = 300, Fm = 200, Fs.n = 30, and n = 10®®cm“®i?“^ 
(i.e. At = 1/3), respectively. The parameters Ce = 0.1, 
eb = 0.01 and p = 2.3 are adopted for the FS and RS 
shocks as well as the trail. Where Etot is the total energy 
of the initial outflow; A is the width of the initial I-ejecta; 
Fn (Fs.n) is the bulk Lorentz factor of fast (slow) neutrons; 
Fm is the initial bulk Lorentz factor of I-ejecta. (c) R-band 
lightcurves, but the IC cooling effect due to the prompt 
7 —rays overlaping with the shocked region and the trail 
has been taken into account. The averaged 7 —ray lumi¬ 
nosity is taken as = lO^^ergs s“^. Other paramters 
and line styles are the same as those in (b). From ref M- 


sumed that the I-ejecta moves slower than the fast 
neutrons). A gap likely forms between the leading 
trail and the I-ejecta. The former forms a distinct trail 
ejecta (T-ejecta) which interacts with the out trail or 
ISM. The latter catches up later and gives rise to a 
rebrigtening signature. Before collision, the radiation 
is dominated by the forward shock emission. Dur¬ 
ing the collision, both the forward shock emission and 
the refreshed shocks (especially the refreshed reverse 
shock) are important. The unshocked trail emission 
is not important in this case. A typical neutron-rich 
ISM-interaction lightcurve is characterized by a slow 
initial rising lightcurve followed by a prominent bump 
signature around tens to hundreds of seconds (Fig0. 

For all the cases, the predicted signatures can be 



'ob, (S) 

FIG. 6 : The early opticcal afterglow lightcurves of a 
neutron-fed long GRB in the ISM interaction case, (a) 
The dynamical evolution of the region shocked by FS as 
a function of time. The dotted, dash-dotted, dashed, and 
solid lines represent x = 0-0, 0.1, 0.5, 1.0 respectively, 
(b) R-band lightcurves. Line styles are the same as in (a). 
Thick lines represent the total early R-band lightcurves, 
while the thin lines are for RS emission only. Only the RS 
emission for y = 0 and the RRS emission for y = 1 cases 
are plotted. The initial parameters are the same as those 
listed in the caption of Figure 0 except that n = lcm“® 
is adopted. From ref m. 


detected by the UVOT on board the Swift observa¬ 
tory. However, most of these signatures (such as the 
plateau and the bump signature) are not exclusively 
for neutron decay. More detailed modeling and case 
study are needed to verify the existence of the neutron 
component. 


IV. DISCUSSION 

Most of the current afterglow observations take 
place hours after the burst trigger. At this stage, 
the observed afterglow emission are powered by the 
forward shock, so that essentially all the initial in¬ 
formation of the ejecta is lost. In order to diagnose 
the ejecta composition, well-monitored early afterglow 
data are needed, since the early afterglow emission, 
especially in the optical band, is believed to be domi¬ 
nated by the reverse shock emission and possibly the 
trail emission. The reverse shock propagates into the 
ejecta, so the emission property, especially the linear 
polarization degree, depends on the magnetization of 
the ejecta. The neutron signature emission lasts only 
tens to hundreds of seconds, which overlaps with the 
very early shocks emission. Only rapid responce to 
the trigger can catch it. 

In this proceedings paper, we have summarized the 
early afterglow signatures of magnetized GRBs and 
neutron-rich GRBs. These signatures are likely de- 
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tectable by the Ultraviolet Optical Telecope (UVOT) 
on board the Swift observatory. Close monitoring of 
early afterglows from 10s to 1000s of seconds, when 
combined with detailed theoretical modeling, could 
be used to diagnose their existence, which may in turn 
help us to reveal the nature of the GRB central engine 
source. 

Currently, there are just two well studied cases, 
i.e. GRB 990123, and GRB 021211, whose early af¬ 
terglow lightcurves are well consistent with reverse 
shock emission from a moderately-magnetized flow 
11 US HI El m. However, the magnetization it¬ 
self does not mean the outflow is initially magnetized 
since in the internal shock phase, random magnetic 
fields are generated, significant part of which can not 
be dissipated effectively. These magnetic fields could 
be retained in the external shock phase to dominate 
the reverse shock synchrotron emission. In this case. 


there is no net polarization expected. Consequently, 
the early polarization detection is necessary to draw 
definitive conclusions on the initial magnetization of 
the ejecta. 
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